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ABSTRACT 

Measurements  of  wind  speed,  wind  direction,  and  the  vertical  component  of 
turbulence,  from  four  different  commercially  available  Doppler  sodars,  are 
compared  with  similar  measurements  from  in  situ  sensors  on  a  300  m  instru- 
mented tower.  Results  indicate  that  the  four  sodars  measure  wind  speed  and 
direction  accurately  and  with  reasonably  high  precision.  The  sodars  tended  to 
overestimate  the  vertical  component  of  turbulence  at  night  and  to  underesti- 
mate it  during  the  day.  Precision  in  those  measurements  was  considerably 
poorer  than  for  the  averaged  speeds  and  directions.  Analysis  of  the  vertical 
wind  from  the  sodars  indicates  that  the  measurement  inaccuracies  arise  from  a 
combination  of  aliasing  and  spatial  averaging. 
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1.   INTRODUCTION 

During  the  first  three  weeks  of  September  1982,  an  experiment  was  con- 
ducted at  the  Boulder  Atmospheric  Observatory  (BAO) ,  under  the  sponsorship  of 
the  Environmental  Protection  Agency  (EPA),  to  assess  the  ability  of  in  situ 
and  remote  sensors  to  measure  the  mean  and  turbulent  properties  of  the  lower 
atmosphere.  The  experiment  was  conducted  in  response  to  the  need  for  com- 
parative data  from  which  scientists  could  evaluate  the  accuracy,  precision, 
and  general  performance  of  some  of  the  more  commonly  used  meteorological 
instruments  that  measure  atmospheric  turbulence. 

Recent  advances  in  modeling  transport  and  diffusion  of  pollutants, 
achieved  largely  through  theoretical  insights  gained  from  field  experiments, 
point  to  the  site-specific  nature  of  turbulence.  Attention  is  therefore  being 
directed  to  better  on-site  characterization  of  turbulence  and  to  the  develop- 
ment of  techniques  for  measuring  the  mean  and  turbulent  wind  variables  needed 
for  input  into  the  models.  This  experiment  was  designed  to  provide  infor- 
mation needed  to  formulate  a  monitoring  strategy  for  developing  site-specific 
dispersion  meteorology. 

The  Boulder  Atmospheric  Observatory  (BAO)  was  chosen  as  the  site  for  the 
experiment  because  of  the  availability  of  precise  profile  and  turbulence  data 
from  sensors  on  a  300  m  tower.  Facilities  for  launching  rawinsondes  and  for 
processing  the  data  received  from  them  were  added  benefits.  Two  categories  of 
sensors  were  tested  against  measurements  on  the  tower.  One  category  consisted 


of  a  set  of  lightweight  in  situ  sensors  selected  for  this  study.  The  other 
category  consisted  of  four  commercially  available  Doppler  sodars  with  the 
capability  to  measure  variances  in  the  vertical  wind  component  in  addition  to 
the  mean  three-dimensional  wind  field.  This  report  deals  only  with  the  sodar 
comparisons.  The  in  situ  sensor  comparisons  will  be  described  in  a  separate 
report.  The  question  addressed  here  is  whether  the  sodars  can  measure  the 
mean  and  turbulent  properties  of  the  flow  at  heights  that  are  above  the  reach 
of  conventional  10  m  meteorological  towers. 


2.  DESCRIPTION  OF  INSTRUMENTATION 

Four  Doppler  sodar  manufacturers  who  currently  market  their  products  in 
the  United  States  were  invited  by  EPA  to  participate  in  the  experiment.  Under 
arrangements  made  through  EPA's  principal  contractor,  Meteorology  Research, 
Inc.,  the  sodars  were  installed  and  operated  by  personnel  from  the  par- 
ticipating firms.  The  four  systems  differ  significantly  in  their  physical 
configuration  and  approach  to  signal  processing.  Two  of  the  systems  used 
monostatic  three-axis  arrays,  one  used  a  bistatic  array,  and  the  other  a  colo- 
cated  monostatic/bi static  array. 

2.1  AeroVironment  Three-Axis  Monostatic  System  (AV) 

The  system  consisted  of  three  acoustic  transceivers  mounted  on  a  trailer. 
Doppler  shifts  in  the  backscattered  signals  received  on  each  axis  were 
interpreted  as  wind  components  in  the  radial  directions.  Wind  components  thus 
measured  were  transformed  into  components  along  the  N-S,  E-W,  and  vertical 
directions.  Since,  in  this  configuration,  sampling  volumes  are  separated  by 
large  distances  the  assumption  of  horizontal  homogeneity  in  the  mean  wind 
field  is  essential  to  justify  using  wind  components  measured  along  the  dif- 
ferent axes  for  the  coordinate  transformation. 

The  AeroVironment  system  transmitted  a  sound  pulse  (150-200  watts)  at  a 
frequency  of  1500  Hz  (duration  0.18  s)  sequentially  from  each  of  three  adja- 


cent  pencil -beam  antennas.  One  tilts  south  30°  from  the  vertical  to  be  sen- 
sitive to  the  N-S  component,  one  tilts  west  30°  from  vertical  to  be  sensitive 
to  the  E-W  component,  and  one  points  straight  up  to  be  sensitive  to  the  ver- 
tical component.  The  receiver  echo  is  heterodyned  and  then  passed  through  an 
electronic  comb  filter  with  31  teeth  to  yield  the  full  spectral  distribution 
in  the  return  signal.  For  each  33.3  m  altitude  range  gate,  the  spectrum  is 
examined  according  to  several  criteria  to  obtain  a  best  estimate  of  Doppler 
shift  along  with  an  estimated  reliability  factor.  The  pulse  repetition  inter- 
val was  8  s. 

2.2  Remtech  Three-Axis  Monostatic  System  (REM) 

The  Remtech  system  (developed  originally  at  Bertin  et  Cie)  also  uses  a 
trailer-mounted  array  of  three  transceivers.  They  are  operated  in  sequence  as 
monostatic  systems;  the  same  assumption  of  horizontal  homogeneity  is  invoked 
for  wind  measurements.  In  this  system  the  horizontal  wind-sensing  antennas 
are  tilted  18°  from  the  vertical;  the  transmitted  pulse  is  a  1600  Hz  signal  of 
0.08  s  duration.  The  received  signals  are  digitized  after  appropriate  band- 
pass filtering,  and  the  Doppler  frequency  shift  extracted  by  using  Fast 
Fourier  Transform  (FFT)  techniques.  The  pulse  repetition  interval  was  5  s. 

2.3  Radian  Corporation  Colocated  Monostatic/Bi static  System  (RAD) 

Radian's  antenna  configuration  permitted  both  monostatic  and  bistatic 
operation.  Both  systems  shared  the  central,  vertically  pointing,  pencil -beam 
transceiver.  The  two  tilted  (18°  from  vertical)  monostatic  transceivers  were 
not  located  close  to  the  vertical  transceiver  as  on  the  AV  and  REM  systems  but 
aimed  to  intersect  at  a  height  of  150  m.  In  the  bistatic  mode,  two  fan-beam 


transmitters  (located  250  m  to  the  south  and  to  the  west)  illuminated  the  ver- 
tical beam  of  the  central  transceiver.  The  movement  of  the  sound  pulse  up  the 
vertical  beam  was  followed  by  time  gating  of  the  receiver  signal.  Ooppler 
frequency  shifts  in  each  gated  segment  are  converted  to  wind  velocity  com- 
ponents to  produce  a  wind  profile.  For  details  of  this  approach  see  Kaimal 
and  Haugen  (1977). 

In  both  configurations,  the  vertical  transceiver  was  operated  in  the 
monostatic  mode  to  measure  the  vertical  wind  component.  Since  the  three 
monostatic  beams  were  not  divergent,  the  assumption  of  homogeneity  is  not  as 
critical  here.  The  RAD  system  transmitted  120  watt  pulses  at  2.0  kHz  (0.1  s 
duration)  and  computed  Doppler  shifts  using  the  Complex  Covariance  method. 
RAD  operated  in  three  modes,  monostatic,  bistatic,  and  multistatic 
(alternating  in  time,  one  series  of  monostatic  and  one  series  of  bistatic 
pulses).  Pulse  repetition  interval  for  all  systems  was  5  s. 

2.4  Xontech  Three-Axis  Bistatic  System  (XON) 

This  bistatic  system  consisted  of  a  vertical  pencil -beam  transceiver  and 
two  fan-beam  receivers  aimed  at  a  central  vertical  common  volume.  Therefore, 
the  geometry  in  this  system  is  exactly  the  reverse  of  the  RAD  bistatic  system. 
The  transmitted  frequency  was  2.0  kHz  (0.08  and  0.16  s  duration  under  computer 
control).  Its  bistatic  baseline  was  350  m  long.  The  fan-beam  antennas 
receive  signals  scattered  from  the  vertical  transceiver  beam,  so  the  winds 
are  computed  along  a  vertical  column  above  the  transceiver  (Balser  et  al., 
1976)  as  in  the  RAD  bistatic  system.  A  microcomputer  determined  the  Doppler 
frequency  shift  with  an  FFT  detection  scheme  powerful  enough  to  sense  small 
frequency  variations  in  the  presence  of  high  ambient  noise  levels.  The  wind 


data  are,  therefore,  presented  without  qualifiers,  but  if  the  program  cannot 
detect  a  consistent  signal  for  the  entire  averaging  period,  no  data  are 
printed  for  that  height  range.  The  pulse  repetition  interval  was  5  s. 

2.5  BAO  Instrumented  Tower 

The  300  m  tower  at  BAO  is  instrumented  at  eight  levels:  10,  22,  50,  100, 
150,  200,  250,  and  300  m.  Sonic  anemometers  installed  at  each  level  measured 
the  three-dimensional  wind  field.  R.M.  Young  propeller-vane  anemometers  are 
mounted  on  the  side  of  the  tower  opposite  the  sonic  anemometers  to  serve  as 
backup  wind  sensors  when  the  tower  is  shadowing  the  sonic  anemometers.  For 
this  experiment,  the  sonic  anemometers  were  mounted  on  the  booms  pointing  SSW, 
and  the  propeller-vane  anemometers  were  on  booms  pointing  NNE.  These  booms 
also  supported  sensors  for  measuring  mean  and  fluctuating  air  temperatures  and 
the  dewpoint  temperature.  Data  from  the  sonic  anemometer  and  other  fast- 
response  sensors  were  sampled  10  times  per  second;  the  propeller-vane  anemom- 
eter, like  other  slow-response  sensors,  was  sampled  only  once  per  second. 

Data  from  the  BAO  are  recorded  in  one  of  two  modes.  In  the  "regular" 
mode,  only  the  10  s  averaged  data  points  and  10  s  grab  samples  (last  point  in 
a  10  s  data  block)  of  the  time  series  are  retained.  In  the  "raw  data"  mode 
all  data  points  are  recorded.  In  both  modes  the  software  computes  and  lists 
once  every  20  min  the  means,  variances,  and  fluxes  for  the  preceding  20  min 
period.  These  listings  became  the  common  reference  for  comparing  the  perfor- 
mance of  the  different  sodars.  The  raw  data  mode  is  employed  only  when  the 
full  time  series  is  needed  for  special  analyses  or  for  the  details  in  the 
structure  of  the  flow.  Such  recordings  were  made  for  three  relatively  brief 
periods  during  this  experiment  for  the  purpose  of  comparing  the  spectral 


responses  of  the  sodars  and  examining  the  limitations  imposed  by  sampling 
volumes  and  sampling  rates. 

2.6  Rawinsonde 

The  rawinsonde  station  located  near  the  base  of  the  tower  was  operated 
daily  during  this  experiment  to  obtain  wind  data  above  300  m.  In  addition  to 
providing  basic  meteorological  information,  the  data  were  needed  to  assess  the 
range  limitations  of  the  monostatic  sodars  under  different  atmospheric  con- 
ditions. Standard  ground  monitoring  equipment  and  radiosonde  packages  of  the 
type  used  by  the  National  Weather  Service  were  deployed.  The  procedure 
followed  was  to  schedule  a  release  at  about  sunrise  and  about  midday,  such 
that  a  release  coincided  with  the  observing  period  of  a  different  sodar  system 
each  day.  The  purpose  was  to  provide  a  near-equal  distribution  of  radiosonde 
comparisons  for  each  type.  For  two  4  h  periods  during  the  experiment,  hourly 
ascents  were  made  covering  sunrise  to  late  morning  and  late  afternoon  to  sun- 
set. 


3.  DESCRIPTION  OF  THE  FIELD  PROGRAM 

The  sodar  measurements  were  centered  over  an  area  0.5  x  0.3  km,  about  0.65 
km  southwest  of  the  BAO  tower.  The  bistatic  arrays  were  laid  out  to  provide 
a  height  range  of  at  least  300  m,  hence  the  requirement  for  such  a  large  test 
area.  The  deployment  of  the  antennas  for  the  four  systems  is  shown  in  the 
plot  plan  of  Fig.  1.  The  electronic  equipment  associated  with  the  Doppler 
systems  was  housed  in  trailers  located  within  the  visitors'  area.  A  larger 
trailer  in  the  same  area  served  as  the  control  center  for  the  experiment. 

The  terrain  in  the  vicinity  of  the  tower,  including  the  area  covered  in 
Fig.  1,  is  reasonably  flat.  There  is  a  gentle  slope  to  the  north,  east,  and 
west  (see  Fig.  2),  but  the  steepest  grade  (7%)  is  toward  the  small  hill 
south  of  the  tower.  Except  for  the  trailers  and  the  fence  surrounding  the 
visitors'  area,  the  site  is  free  of  small-scale  surface  obstructions.  The 
photographs  in  Figs.  3(a)  and  (b)  show  the  characteristics  of  the  test  site. 
For  a  more  detailed  description  of  the  site  and  its  effect  on  flow  over  the 
area  see  Kaimal  et  al .  (1982). 

Procedures  for  data  collecting  and  reporting  were  established  to  ensure 
against  unfair  bias  for  any  of  the  participants.  All  systems  were  assumed 
capable  of  unattended  continuous  operation.  All  systems  were  to  provide  data 
in  the  form  of  wind  speeds,  wind  directions,  vertical  wind  components,  and 
standard  deviations  of  vertical  wind  averaged  over  20  min  periods  coincident 
with  the  BAO  averaging  periods.  The  three  comparison  levels  would  be  100, 
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Figure  1.  Plot  plan  showing  sodar  antenna  deployment  in  the  BAO  test 
area. 
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Figure  2.  Contour  map  of  the  immediate  BAO  terrain  showing  location  of 
the  tower  and  the  sodar  test  area. 


11 


n*»  WW  ""'■■pj1' 


Figure  3.  (a)  View  of  the  sodar  test  area  looking  north  from  County 
Road  8. 


(b)  Instrumentation  on  the  BAO  300  m  tower, 
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200,  and  300  m.  Every  morning  at  0800  MST  the  data  collected  over  the  pre- 
vious 24  h  period  were  submitted  to  EPA  personnel  directing  the  experiment  in 
exchange  for  tower  data  covering  the  same  period. 

Concurrent  operation  of  some  of  the  sodars  was  considered  at  one  time,  but 
quickly  ruled  out  because  of  cross-contamination,  even  between  systems 
operating  at  different  frequencies.  The  sodars  were  therefore  operated  in 
sequence,  the  switch-over  from  one  system  to  the  other  being  controlled  by  a 
central  timer  switch.  The  assigned  observing  period  was  one  20  min  period 
each  cycle.  The  experiment  took  place  between  1  and  21  September.  On  three 
different  occasions,  the  observing  period  was  extended  to  120  min  (9  Septem- 
ber) and  to  80  minutes  (16/17  and  18/19  September)  in  order  to  obtain  long 
enough  records  for  spectrum  analysis.  The  BA0  data  were  recorded  in  the  raw 
data  mode  on  these  occasions.  The  three  occasions  were:  0800-1600,  9 
September;  1520-0800,  16/17  September  and  1600-0140,  18/19  September.  The  two 
4  h  periods,  when  rawinsonde  measurements  were  made  hourly,  occurred  on  8 
September  (0400-0800)  and  18  September  (1200-1600). 

AV,  REM,  and  RAD  computed  o ,  the  standard  deviation  of  w  (the  vertical 
wind  component),  from  their  time  series.  Missing  data  points  were  not  filled 
in  by  interpolation,  but  the  number  of  points  missed  (or  accepted)  was 
displayed.  REM  used  four-point  block  averages  instead  of  the  original  time 
series.  X0N  computed  its  standard  deviation  from  the  width  of  its  2  min  w 
spectra,  estimated  for  each  level.  Successive  2  min  standard  deviations  were 
averaged  to  obtain  the  20  min  values.  Each  spectrum  was  automatically  exa- 
mined for  level  and  shape  of  background  noise,  and  steps  were  taken  to  remove 
their  effects. 
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AV,  REM,  and  RAD  reported  standard  deviations  of  the  wind  direction,  aD. 
REM  and  RAD  calculated  theirs  from  the  wind  direction  time  series;  AV  used  the 
relationship 

aD  =  av/U  , 

where  U  is  mean  horizontal  wind  vector.  XON  reported  the  standard  devia- 
tions of  the  longitudinal  and  lateral  wind  components.  These  were  obtained 
through  coordinate  transformation  of  horizontal  wind  components  measured  by 
the  sodar. 

No  attempt  is  made  in  this  report  to  present  the  results  of  our  an  com- 
parisons. The  azimuth  direction  standard  deviations  show  very  large  scatter. 
The  data  are  withheld  pending  a  better  understanding  of  the  reasons  for  the 
scatter.  Meanwhile,  we  can  only  suggest  caution  in  using  on  for  diffusion 
predictions. 

During  the  experiment,  the  sodars  encountered  a  wide  range  of  weather 
conditions:  from  clear  skies  to  heavy  rain,  and  winds  ranging  from  very  light 
to  well  over  10  m/s.  A  summary  of  daily  weather  conditions  for  the  duration 
of  the  experiment  is  given  in  the  Appendix.  As  the  manufacturers  were  asked 
to  submit  only  data  they  believed  were  correct,  all  submitted  data  were  used 
in  our  comparison  study. 

AV,  REM,  and  XON  maintained  a  consistent  operating  pattern  throughout  the 
experiment.  However,  RAD  changed  its  operating  mode  every  24  h,  switching 
from  multimode  to  bistatic  and  monostatic,  back  to  multimode  and  so  on.  For 
the  extended  observing  periods,  RAD  operated  in  the  monostatic  mode  on  9 
September,  in  the  bistatic  mode  on  16/17  September,  and  multimode  on  18/19 
September. 
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4.  MEASUREMENT  OF  THE  STANDARD  DEVIATION  OF  w 

At  the  three  heights  under  consideration  (100  m,  200  m,  300  m),  only  AV 
attained  a  complete  data  record.  The  other  three  manufacturers  had  these 
ranges  of  completeness  at  the  three  heights:  RAD,  70%  to  95%;  REM,  55%  to 
72%;  XON,  77%  to  94%.  However,  wind  shadow  effects  depleted  the  sonic  records 
so  that  the  final  outcome  regarding  completeness  shows  composite  sonic/sodar 
percentages  of  47%  at  100  m,  41%  at  200  m,  and  47%  at  300  m  for  the  sodar  ver- 
sus sonic  comparison. 

4.1  Sodar  Reference  Differences 

Since  the  sonic  measurements  of  vertical  wind  speed  standard  deviation 
provide  reference  values,  the  accuracy  and  precision  of  each  sodar  system  can 
be  determined  from  the  collections  of  20  mi n  average  differences.  The  two 
input  variables  for  these  computations  are  a  for  the  sodar  instruments  and 


/ww  for  the  sonic  vertical  wind  speed.  The  comparative  statistics  used  to 
estimate  accuracy  and  precision  then  become  the  average  difference,  or  sample 
bias  (b),  and  the  standard  deviation  of  the  differences  (s).  In  addition,  the 
root  mean  square  difference,  or  comparability  (c),  is  computed;  this  statistic 
was  defined  by  Hoehne  (1971,  1977),  and  it  characterizes  the  repeatability  of 
a  system.  Finally,  the  precision  is  also  represented  as  a  percentage  (s1)  of 
the  average  value  of  sonic  a i.e.,  a  coefficient  of  variation. 
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Expressions  for  b,  c,  s,  and  s'  are  given  in  Eqs.  (1)  to  (4). 

b  =  £  I  [(oW)1  -  Mwn.]  (1) 

c   =   {J-Ji  [(^  -  v^.]2}172  (2) 

s  =  (c2  -  b2)1/2  (3) 

S'  =  s/aw  x  100.  (4) 

Values  of  b,  c,  s,  and  s'  are  presented  in  Table  1  for  the  combined  sodar 
observations  at  each  of  three  heights,  as  well  as  for  the  individual  vendor 
data  subsets. 

The  sample  bias  (b)  in  Table  1  shows  a  large  range  of  values  around  nearly 
constant  composite  values  of  0.08  to  0.09  m/s.  At  100  m  the  spread  is 
greatest,  REM  having  the  only  negative  bias  value  (i.e.,  sodar  <  sonic)  and 
X0N  having  a  sizable  0.23  m/s  bias.  AV  is  well  below  the  composite  bias 
amount;  RAD  is  above  it.  At  200  m  the  REM  value  is  slightly  negative,  AV 
remains  small,  X0N  approximates  the  composite  value,  but  RAD  is  in  excess  of 
0.2  m/s.  At  300  m  the  RAD  bias  continues  to  be  relatively  large,  but  the 
other  vendors  are  grouped  between  0.04  and  0.07  m/s. 

From  Table  1  it  is  clear  that  there  is  much  scatter  in  c  and  s'  about  the 
true  value  in  all  systems.  There  is  no  statistical  difference  between  s  and 
s'  values  for  AV  and  REM  at  any  of  the  levels.  The  division  of  Table  1  into 
daytime  and  nighttime  categories  is  displayed  in  Table  2.  The  bias  columns 
show  that  all  sodar  systems  tend  to  overestimate  a  at  night. 
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Table  1.  Sodar  a  compared  with  sonic  a, 
w  w 


Height     Vendor     b  (m/s)     c  (m/s)     s  (m/s)     s'(%)    N 


100  m 


200  m 


300  m 


All 

0.08 

0.24 

0.22 

SO 

678 

RAD 

0.12 

0.25 

0.21 

47 

178 

REM 

-0.05 

0.18 

0.17 

38 

139 

AV 

0.01 

0.16 

0.16 

35 

190 

XON 

0.23 

0.34 

0.24 

53 

171 

All 

0.08 

0.27 

0.26 

54 

576 

RAD 

0.22 

0.39 

0.32 

65 

144 

REM 

0.00 

0.19 

0.19 

39 

119 

AV 

0.03 

0.20 

0.20 

43 

167 

XON 

0.08 

0.25 

0.24 

51 

146 

All 

0.09 

0.27 

0.26 

54 

665 

RAD 

0.23 

0.38 

0.30 

62 

158 

REM 

0.07 

0.23 

0.22 

47 

136 

AV 

0.04 

0.25 

0.25 

53 

214 

XON 

0.04 

0.19 

0.18 

38 

157 

b  =  bias  (accuracy)  (sodar-sonic) 

c  =  comparability 

s  =  standard  deviation  of  differences  (precision) 

s'  =  s  expressed  as  a  percentage  of  average  value  of  sonic  standard  deviation 

N  =  number  of  observations 


AV  =  Aerovi  ronment 
RAD  =  Radian 
REM  =  Remtech 
XON  =  Xontech 
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Table  2.     Separation  of  Table  1   into  daytime  and  nighttime  categories 


Hei  ght 

Vendor 

bd 

bn 

Cd 

c 
n 

Sd 

s 

n 

S'd 

s' 

n 

Nd 

N 
n 

100  m 

All 

0.02 

0.17 

0.21 

0.27 

0.21 

0.21 

34 

99 

389 

288 

RAD 

0.08 

0.18 

0.23 

0.26 

0.22 

0.20 

34 

88 

100 

78 

REM 

-0.10 

0.22 

0.21 

0.24 

0.19 

0.11 

28 

58 

77 

62 

AV 

-0.03 

0.08 

0.16 

0.15 

0.16 

0.13 

27 

57 

115 

75 

XON 

0.13 

0.37 

0.25 

0.43 

0.21 

0.22 

33 

99 

97 

73 

200  m 

All 

0.04 

0.14 

0.30 

0.23 

0.30 

0.19 

43 

78 

306 

269 

RAD 

0.21 

0.23 

0.43 

0.33 

0.38 

0.24 

54 

93 

78 

66 

REM 

-0.06 

0.06 

0.21 

0.17 

0.20 

0.16 

27 

72 

62 

57 

AV 

-0.02 

0.09 

0.23 

0.18 

0.23 

0.15 

34 

62 

90 

77 

XON 

0.00 

0.17 

0.28 

0.23 

0.28 

0.15 

42 

62 

76 

69 

300  m 

All 

0.06 

0.12 

0.28 

0.27 

0.27 

0.24 

41 

89 

359 

304 

RAD 

0.21 

0.24 

0.39 

0.36 

0.33 

0.27 

49 

98 

87 

71 

REM 

0.03 

0.11 

0.25 

0.22 

0.25 

0.19 

36 

78 

70 

66 

AV 

-0.00 

0.09 

0.23 

0.29 

0.23 

0.27 

36 

100 

121 

93 

XON 

0.03 

0.06 

0.22 

0.15 

0.22 

0.14 

33 

52 

81 

74 

Subscript  d: 
Subscript  n: 


daytime   (0600-1800  hours) 

nighttime   (0000-0600  and   1800-2400  hours) 
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4.2  Individual  20  Minute  Average  Values 

Additional  information  about  the  effectiveness  of  sodar  measurements  of 

o  can  be  sought  in  the  scatter  plots  of  their  20  min  average  values  against 
w 

the  sonic  standard  deviation  of  vertical  wind  speed.  Such  plots  are  presented 
by  height  and  vendor  in  Figs.  4-15;  daytime  observations  (0600-1800  MST)  are 
distinguished  from  nighttime  observations. 

Each  chart  has  a  broken  line  at  45°  from  the  origin  to  represent  a  slope 
of  1.  The  estimates  of  the  correlation  coefficient  based  upon  the  20  day 
sample  are  given  in  Table  3. 

Each  chart  also  has  a  line  for  sodar  regressed  linearly  upon  reference 
according  to  the  following  model: 


where 


Yi  ■  b0  +  e1x1  ♦  «,  , 


Y-  =  ith  sodar  measurement, 

X.  =  ith  reference  measurement, 

8  =  Y-intercept, 
o 

81  =  slope, 

e.  =  error  term, 

In  all  cases,  the  regression  line  misses  the  origin  on  the  upper  side  and 
has  a  slope  less  than  1.  Both  characteristics  are  significant  at  the  0.05 
level . 

As  is  readily  seen  in  Figs.  4-15,  the  sodar  systems  are  all  consistently 
recording  too  high  when  aw  <  0.5  m/s  and  too  low  when  a  >  1.0  m/s. 
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Table  3.  Regression  analysis  for  «. 

w 


Height        Vendor 


100  m 


200  m 


300  m 


RAD 

0.85 

0.28 

0.67 

178 

REM 

0.89 

0.10 

0.68 

139 

AV 

0.89 

0.18 

0.63 

190 

XON 

0.67 

0.46 

0.51 

171 

RAD 

0.77 

0.39 

0.67 

144 

REM 

0.88 

0.12 

0.75 

119 

AV 

0.84 

0.20 

0.64 

167 

XON 

0.83 

0.26 

0.63 

146 

RAD 

0.72 

0.39 

0.69 

158 

REM 

0.84 

0.17 

0.79 

136 

AV 

0.70 

0.19 

0.68 

214 

XON 

0.82 

0.13 

0.80 

157 

p     =  estimate  of  correlation  coefficient 

3q  =  intercept  term 

3,   =  slope  term 

N     =  number  of  observations 
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4.3  Sodar  Modes 

The  sodar  systems  had  two  distinct  operational  modes:  monostatic  (MONO) 
and  bistatic  (BI).  AV  and  REM  employed  the  former  mode;  XON  represented  the 
latter  mode;  RAD  rotated  daily  among  MONO,  BI,  and  a  combination  of  the  two 
modes  called  multistatic  (MULTI). 

A  comparison  of  RAD  modal  data  with  simultaneous  sonic  standard  deviation 
values  yields  the  b,  c,  and  s  results  of  Table  4.  The  sample  bias  (b)  is 
significantly  nonzero  in  all  modes  at  all  heights.  It  is  equivalent  among 
modes  at  each  height  at  a  probability  level  of  0.05,  except  for  a  significant 
difference  between  MONO  and  MULTI  biases  at  200  m. 

Comparability  (c)  and  standard  deviation  (s)  do  not  show  a  uniform  ranking 

of  modes  with  height.  The  large  c-value  for  MONO  at  200  m  is  partly  due  to 

the  large  bias  and  the  presence  of  two  a     values  slightly  greater  than  2  m/s 

in  the  subset.  The  s-values  are  inconsistent  in  their  equivalence  from  height 

to  height,  possibly  because  of  the  sparseness  of  the  data. 

Scatter  diagrams  of  individual  20  min  average  values  of  sodar  a  versus 
J  3  w 

sonic  a  are  shown  by  height  and  by  mode  in  Figs.  16-24.  The  BI  mode  has  the 
w 

largest  departures  from  the  1:1  line.  Figures  16-18  do  not  include  sonic  a  > 

1.0  m/s.  Thus  the  plots  involving  BI  mode  are  not  strictly  comparable  with 

the  others. 

In  summary,  there  is  a  consistent  tendency  in  all  sodar  systems  to 

overestimate  o    at  low  sonic  readings  (in  stable  conditions).  Conversely,  the 
w 

sodar  systems  underestimated  o    when  the  true  value  was  higher  (unstable 
conditions).  In  all  cases  the  regression  lines  had  positive  y-intercepts  and 
slopes  less  than  unity. 
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On  average,  sample  biases  were  generally  significantly  positive,  but  AV 
did  not  show  a  significant  bias  at  100  m  and  was  marginally  biased  at  200  m 
and  300  m.  REM  was  the  only  vendor  with  a  negative  bias  at  100  m,  but  at  300 
m,  REM  had  a  positive  bias,  as  did  all  the  others. 

REM  performed  very  well  with  respect  to  c  and  s  at  200  m,  but  results  for 
these  statistics  were  not  consistent  with  height.  In  fact,  a  different  vendor 
emerged  with  the  lowest,  i.e.,  best,  values  of  c  and  s  at  each  level. 

Mode  switching  in  the  RAD  system  revealed  a  significant  amount  of  bias  in 
each  mode  (BI,  MONO,  MULTI)  at  100  m,  200  m,  and  300  m.  However,  as  Table  4 
indicates,  there  was  little  difference  among  modes. 
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Table  4.  Radian  modal  a     compared  with  sonic  a. 

w  w 


Height         Mode         b  (m/s)     c  (m/s)       s  (m/s) 


100  m 


200  m 


300  m 


MONO 

0.14 

0.28 

0.24 

68 

BI 

0.14 

0.24 

0.19 

44 

MULTI 

0.10 

0.22 

0.20 

66 

MONO 

0.30 

0.50 

0.40 

58 

BI 

0.19 

0.35 

0.29 

30 

MULTI 

0.15 

0.26 

0.21 

56 

MONO 

0.25 

0.36 

0.25 

57 

BI 

0.23 

0.40 

0.32 

48 

MULTI 

0.19 

0.39 

0.34 

53 

b  =  bias  (accuracy) 

c  =  comparability 

s  =  standard  deviation  of  differences  (precision) 

N  =  number  of  observations 
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5.  MEASUREMENT  OF  WIND  SPEED 

Measurements  of  wind  speed  (S)  were  obtained  at  sampling  rates  of  1 
datum/15  s  (REM,  RAD,  XON)  and  1  datum/24  s  (AV),  and  average  S  was  computed 
in  the  field  over  20  min  intervals.  The  four  sodars  cycled  in  sequence,  so 
that  a  maximum  of  about  1440/4,  or  360,  values  could  be  obtained  by  each  ven- 
dor in  the  20  days  of  the  experiment.  At  the  three  heights  under  con- 
sideration (100  m,  200  m,  300  m),  AV  had  complete  data  and  the  other 
manufacturers  ranged  in  completeness  as  follows:  RAD,  60%  to  96%,  REM,  55%  to 
72%;  XON,  91%  to  92%,  depending  on  the  variable. 

To  examine  the  accuracy  and  precision  of  the  sodars,  simultaneous  obser- 
vations of  wind  speed  were  recorded  from  sonic  anemometers  at  the  same  three 
heights  on  a  tower  that  was  about  600  m  from  the  sodar  systems.  The  sonic 
systems  had  a  sampling  rate  of  10  Hz,  and  they  are  regarded  as  the  reference 
instruments  in  the  evaluation.  However,  owing  to  a  wind  shadow  zone  created 
by  the  tower,  extending  ±40°  from  north  for  the  sonic  instruments,  reference 
data  in  this  sector  were  obtained  by  the  propeller-vane  at  the  BA0  tower.  A 
comparison  of  sonic  and  propeller  wind  speed  measurements  on  the  tower  showed 
that  the  instruments  were  approximately  equivalent.  The  resulting  sonic  and 
propeller  data  sets  are  about  equal  in  size  at  200  m,  but  the  sonic  set  has 
one-third  more  data  at  100  m  and  four-thirds  more  data  at  300  m.  Considering 
omissions  in  the  reference  data,  these  completeness  percentages  resulted:  AV, 
83%  to  91%;  RAD,  55%  to  88%;  REM,  51%  to  66%;  XON,  75%  to  84%. 
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5.1  Sodar  Reference  Differences 

Values  of  sample  bias  (b),  comparability  (c),  and  standard  deviation  (s) 
and  coefficient  of  variation  (s1)  for  the  differences  between  sodar  and 
reference  values  are  presented  in  Table  5  for  combined  sodar  observations  at 
each  height  as  well  as  for  the  sodar  record  of  each  vendor.  Propeller  wind 
speeds  were  excluded  when  the  wind  speed  was  less  than  1  m/s. 

The  estimates  of  bias  in  Table  5  show  mostly  negative  values  at  100  m 
and  a  composite  value  near  -0.4  m/s.  Since  the  difference  is  taken  as  (sodar 
-  reference),  this  means  that  the  sodar  systems  tend  to  register  too  low.  An 
exception  is  RAD,  which  does  not  have  a  significant  bias  at  100  m. 

At  200  m,  the  vendors  all  record  too  high,  and  at  300  m,  RAD  and  X0N  again 
record  too  high  whereas  AV  is  slightly  negative  and  REM  unbiased.  Biases  were 
computed  for  day  (0600-1800  hours)  and  night  (1800-0600  hours)  values.  Most 
differences  between  day  and  night  are  insignificant. 

The  comparability  (c)  of  sodar  wind  speeds  with  reference  values  is  also 
given  in  Table  5.  Precision  is  represented  by  standard  deviation  (s)  and  per- 
centage deviation  (s1).  The  s'  values  range  from  about  15%  to  35%  around  com- 
posite values  near  25%. 

5.2  Individual  20  Minute  Average  Values 

Additional  information  about  the  characteristics  of  sodar  measurements  of 
S  can  be  sought  in  the  scatter  diagrams  of  sodar  20  mi n  average  values  plotted 
against  reference  values.  Such  plots  are  presented  by  height  and  by  vendor  in 
Figs.  25-36. 

Each  chart  has  a  broken  line  at  45°  from  the  origin  representing  a  slope 
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of  1.  The  estimates  of  the  correlation  coefficient,  slope,  and  intercept  are 
given  in  Table  6. 

The  agreement  between  sodar  and  tower  wind  speed  measurements  is  obviously 
quite  good.  Differences  between  manufacturers  can  be  deduced  by  the  reader. 

5.3  Sodar  Modes 

A  comparison  between  RAD  modal  data  and  corresponding  20  min  reference 
averages  yields  the  b,  c,  and  s  results  of  Table  7.  The  monostatic  mode  had 
significantly  higher  bias  than  the  other  two  modes  at  all  heights. 

The  comparability  and  standard  deviation  show  a  distinct  advantage  to  the 
bistatic  mode;  MONO  has  the  greatest  magnitudes  of  c  at  200  m  and  300  m,  and 
MULTI  has  the  greatest  magnitudes  of  s  at  100  m  and  300  m. 
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Table  5.  Sodar  wind  speed  compared  with  reference  wind  speed 


Height     Vendor     b  (m/s)     c  (m/s)     s  (m/s)     s'(%)    N 


100  m 


200  m 


300  m 


All 

-0.42 

1.28 

1.21 

28 

1179 

AV 

-0.50 

1.03 

0.90 

21 

327 

RAD 

0.02 

1.18 

1.18 

28 

315 

REM 

-0.12 

0.62 

0.60 

14 

236 

XON 

-1.04 

1.88 

1.56 

37 

301 

All 

0.14 

0.98 

0.96 

23 

1019 

AV 

0.05 

0.72 

0.72 

17 

298 

RAD 

0.31 

1.00 

1.47 

35 

258 

REM 

0.12 

0.73 

0.72 

17 

194 

XON 

0.09 

0.71 

0.70 

17 

269 

All 

0.16 

1.24 

1.23 

27 

1005 

AV 

-0.10 

1.15 

1.15 

25 

328 

RAD 

0.29 

1.71 

1.69 

37 

198 

REM 

0.02 

0.74 

0.74 

17 

183 

XON 

0.44 

1.20 

1.12 

25 

296 

b  =  bias  (accuracy) 

c  =  comparability 

s  =  standard  deviation  of  differences  (precision) 

s'  =  s  expressed  as  a  percentage  of  average  value  of  reference  wind  speed 

N  =  number  of  observations 


AV  =  Aerovi  ronment 
RAD  =  Radian 
REM  =  Remtech 
XON  =  Xontech 
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Table  6.  Regression  analysis  for  wind  speed 


Height 


Vendor 


100  m 


200  m 


300  m 


AV 

0.94 

-0.03 

0.89 

327 

RAD 

0.90 

0.41 

0.91 

315 

REM 

0.97 

-0.22 

1.02 

236 

XON 

0.82 

-0.07 

0.77 

301 

AV 

0.96 

-0.02 

1.02 

298 

RAD 

0.87 

0.04 

1.07 

258 

REM 

0.96 

0.34 

0.95 

194 

XON 

0.96 

0.02 

1.01 

269 

AV 

0.93 

0.45 

0.88 

328 

RAD 

0.85 

1.02 

0.84 

198 

REM 

0.96 

0.18 

0.96 

183 

XON 

0.93 

0.46 

0.97 

296 

p  = 
*0  = 

h  = 

N  = 


estimate  of  correlation  coefficient 

intercept  term 

slope  term 

number  of  observations 


51 


Table  7.  Sodar  Radian  modes:  Accuracy  and  precision  for  wind  speed 


Height 


Mode 


b  (m/s) 


c  (m/s) 


s  (m/s) 


100  m 


200  m 


300  m 


BI 

MULTI 

MONO 

-0.15 

-0.10 

0.31 

0.56 
1.54 
1.18 

0.54 
1.53 
1.04 

90 

128 

97 

BI 

MULTI 

MONO 

-0.01 
0.27 
0.74 

0.90 
1.64 
1.79 

0.90 
1.62 
1.63 

78 

110 

70 

BI 

MULTI 

MONO 

0.08 
0.12 
1.01 

1.26 
1.89 
1.96 

1.26 
1.88 
1.68 

66 
92 

40 

b  =  bias  (accuracy) 

c  =  comparability 

s  =  standard  deviation  of  differences  (precision) 

N  =  number  of  observations 
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6.  MEASUREMENT  OF  WIND  DIRECTION 

All  four  sodar  vendors  reported  wind  direction  (D).  At  the  three  heights 
under  consideration  (100  m,  200  m,  300  m),  AV  had  complete  data,  and  the  other 
manufacturers  ranged  in  completeness  as  follows:  RAD,  60%  to  96%;  REM,  55%  to 
72%;  XON,  91%  to  92%.  However,  the  wind  shadow  zone  mentioned  in  Sec.  4 
limited  the  amount  of  data  that  could  be  used.  An  investigation  of  the 
propeller-vane  data  indicated  that  they  could  not  be  substituted  as  reference 
values  for  wind  direction.  There  were  unexplained  disparities  between  the 
propeller  and  sonic  data  that  are  still  under  investigation.  Because  sonic 
data  showed  more  consistent  behavior  and  we  therefore  believe  them  to  be  more 
reliable,  only  the  sonic  wind  direction  measurements  are  used  as  reference 
data,  resulting  in  these  completeness  percentages  for  sodar/sonic  differences: 
AV,  43%  to  63%;  RAD,  36%  to  49%;  REM,  31%  to  39%;  XON,  36%  to  55%. 

6.1  Sodar  Reference  Differences 

Values  of  sample  bias  (b),  comparability  (c),  and  standard  deviation  (s) 
for  the  differences  between  sodar  and  sonic  reference  values  are  presented  in 
Table  8  for  combined  sodar  observations  at  each  height  as  well  as  for  the 
sodar  record  of  each  vendor. 

The  values  of  bias  in  Table  8  show  negative  values  at  100  m  and  200  m  for 
all  vendors,  but  positive  values  at  300  m  for  all  vendors  except  for  RAD. 
However,  most  of  these  values  are  not  significantly  different  from  zero,  con- 


65 


sidering  the  variability  of  the  wind  direction  and  the  number  of  cases 
included.  Day  and  night  biases  given  in  Table  9  show  some  differences,  but 
for  the  most  part  they  are  not  statistically  significant. 

The  comparability  (c)  of  sodar  wind  direction  with  sonic  reference  values 
is  also  given  in  Table  8.  The  low  values  at  200  m  may  be  due  to  the  loss  of 
data  at  that  height  during  a  storm  on  13  September.  Considering  the  scatter 
in  data,  it  can  be  assumed  that  the  vendors'  measurements  of  wind  directions 
are  equivalent  to  each  other. 

6.2  Individual  20  Minute  Averaged  Values 

Scatter  diagrams  of  20  min  averaged  sodar  values  plotted  against  reference 
values  are  given  in  Figs.  37-48.  For  the  most  part,  each  sodar  datum  agreed 
with  the  tower  datum  quite  well.  There  were  some  notable  exceptions,  however, 
especially  during  the  night.  XON  seems  to  have  predicted  many  more  north 
winds  than  actually  occurred  at  100  m  and  300  m,  but  its  agreement  at  200  m  is 
very  good. 

From  the  figures  it  seems  that  most  of  the  vendors  predicted  more 
northerly  winds  at  night  at  100  m  than  actually  occurred,  especially  when  the 
tower  measured  winds  from  the  east  and  south.  This  does  not  seem  to  be  the 
case  at  200  m  or  300  m  and  may  be  due  to  local  terrain  effects,  rather  than 
instrument  problems. 

Correlation  coefficients  are  given  in  Table  10,  with  intercepts  and  slopes 
of  the  linear  regression  lines.  With  two  exceptions,  p  values  are  0.9  or 
greater. 
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6.3  Sodar  Modes 

The  three  modes  of  operation  for  the  RAD  system  showed  no  appreciable  dif- 
ference in  their  ability  to  measure  wind  direction. 


Table  8.  Sodar  wind  direction  compared  with  sonic  wind  direction 


Height       Vendor      b  (deg)       c  (deg)       s  (deg) 


100  m 


200  m 


300  m 


All 

-4.41 

28.59 

28.25 

667 

AV 

-3.87 

26.70 

26.41 

187 

RAD 

-6.76 

27.06 

26.20 

177 

REM 

-2.03 

18.51 

18.40 

137 

X0N 

-4.49 

37.85 

37.58 

166 

All 

-3.43 

23.22 

22.97 

523 

AV 

-0.79 

19.47 

19.45 

155 

RAD 

-7.86 

25.67 

24.43 

128 

REM 

-3.89 

24.67 

24.36 

110 

X0N 

-1.85 

23.80 

23.73 

130 

All 

0.75 

29.59 

29.58 

697 

AV 

0.26 

28.56 

28.56 

227 

RAD 

-3.25 

29.98 

29.80 

131 

REM 

0.62 

19.70 

19.69 

142 

XON 

4.05 

35.96 

35.73 

197 

b  =  bias  (accuracy) 

c  =  comparability 

s  =  standard  deviation  of  differences  (precision) 

N  =  number  of  observations 
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Table  9.  Test  of  day  versus  night  wind  direction  bias 


Height 


Vendor 


b  (day) 


b  (night) 


N  (day) 


N  (night) 


ino  m 


200  m 


300  m 


All 

-6.23 

-2.11 

373 

AV 

-6.81 

0.53 

112 

RAD 

-8.24 

-5.20 

91 

REM 

-3.13 

-0.66 

76 

XON 

-6.09 

-2.41 

94 

All 

-5.75 

-0.84 

94 

AV 

-1.48 

0.06 

86 

RAD 

-11.31 

-4.19 

66 

REM 

-5.25 

-2.43 

57 

XON 

-6.18 

2.89 

68 

All 

-0.44 

2.24 

388 

AV 

-1.01 

2.00 

131 

RAD 

-2.42 

-4.20 

70 

REM 

-0.13 

1.49 

76 

XON 

1.26 

7.66 

111 

294 
75 
86 
61 
72 

72 
69 
62 
53 
62 

309 
96 
61 
66 
86 


Table  10.  Regression  analysis  for  wind  direction 


Height 

Vendor 

p 

eo 

h 

N 

100  m 

AV 

0.94 

-0.25 

0.98 

187 

RAD 

0.92 

5.20 

0.94 

177 

REM 

0.96 

4.26 

0.98 

137 

XON 

0.85 

8.29 

0.96 

166 

200  m 

AV 

0.97 

8.15 

0.95 

155 

RAD 

0.90 

14.07 

0.90 

128 

REM 

0.93 

16.72 

0.91 

110 

XON 

0.95 

2.37 

0.98 

130 

300  m 

AV 

0.93 

-5.86 

1.05 

227 

RAD 

0.93 

-8.87 

1.03 

131 

REM 

0.97 

-4.64 

1.03 

142 

XON 

0.89 

-8.79 

1.09 

197 

p  - 
e0  = 

h  ■ 

N  = 


estimate  of  correlation  coefficient 

intercept  term 

slope  term 

number  of  observations 
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7.  SODAR  RAWINSONDE  COMPARISONS 

Although  our  formal  evaluation  of  the  sodars  was  limited  to  the  three 
observing  levels  (100,  200  and  300  m)  within  the  range  of  the  tower  instru- 
mentation, the  vendors  submitted  data  from  higher  levels  in  the  atmosphere  on 
some  occasions.  These  data  were  compared  with  the  rawinsonde  data  obtained 
concurrently.  Values  of  sample  bias,  comparability,  and  correlation  for  the 
measured  wind  speeds  and  wind  directions  for  each  sodar  are  given  in  Table  11. 
Only  measurements  at  200  m  and  above  are  included  in  the  statistics  because  of 
the  limitations  in  the  rawinsonde's  accuracy  below  that  height.  These  com- 
parisons serve  two  objectives:  to  assess  at  least  in  a  qualitative  sense  the 
ability  of  sodars  to  track  the  broad  features  in  the  speed  and  direction  pro- 
files above  300  m;  and  to  determine  the  magnitude  and  nature  of  the  uncertain- 
ties encountered  in  comparing  measurements  from  a  single  rawinsonde  traverse 
with  time-averaged  (20  min)  sodar  measurements. 

The  numbers  in  Table  11  are  of  the  same  order  as  the  numbers  in  earlier 
tables  comparing  sodar  and  tower  measurements.  Agreement  with  the  rawinsonde 
is  better  for  some  sodars  than  for  others.  The  differences  are  not  considered 
significant,  given  the  fact  that  this  experiment  was  designed,  and  the  anten- 
nas set  up,  for  performance  evaluation  in  the  first  300  m. 

Some  idea  of  the  agreement  in  the  data  can  be  obtained  from  the  scatter 
diagram  in  Fig.  49.  The  wind  speed  and  direction  measurements  from  one  of  the 
sodars  (AV)  is  plotted  against  the  rawinsonde  measurements.  The  scatter  is 


81 


Table  11.  Sodar  wind  speeds  and  directions  for  all  heights 
compared  with  rawinsonde  speeds  and  directions 


Variable 

Sensor 

b 

c 

P 

N 

Speed 

AV 

0.37 

1.17 

0.86 

61 

(m/s) 

RAD 

-0.09 

0.92 

0.61 

14 

REM 

0.24 

0.91 

0.94 

19 

XON 

-1.25 

2.34 

0.38 

31 

Di  recti  on 

AV 

4 

44 

0.90 

61 

(deg) 

RAD 

18 

40 

0.89 

14 

REM 

-  8 

37 

0.95 

19 

XON 

15 

39 

0.96 

31 

b  =  bias 

c  =  comparability 

p  =  estimate  of  correlation  coefficient 

N  =  number  of  observations 
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Figure  49.  Comparison  of  wind  speeds  and  wind  directions  from  rawin^ 
sonde  and  sodar  measurements. 
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not  significantly  larger  than  in  the  plots  against  the  tower  measurements.  It 
is  apparent,  however,  that  the  wind  speeds  agree  better  when  they  are  less 
than  4  m/s,  although  the  percentage  error  is  approximately  the  same.  The 
increase  in  scatter  above  4  m/s  suggests  that  the  larger  spatial  separation 
introduced  between  sensing  volumes  at  the  larger  wind  speeds  (with  the  rawin- 
sonde  drifting  farther  away  from  the  release  point)  is  a  factor  to  be 
recognized  when  rawinsondes  are  used  for  evaluating  sodar  performance. 
Another  factor  to  be  recognized  (but  not  obvious  in  Fig.  49)  is  the  possibi- 
lity of  large  wind  direction  differences  in  sodar  rawinsonde  comparison  under 
light  wind  conditions.  Wind  directions  under  such  conditions  tend  to  be 
highly  variable  both  spatially  and  temporally.  These  conditions  occur  fre- 
quently at  the  BAO  when  the  winds  are  from  E  to  SE.  One  cannot  expect  good 
agreement  between  the  near-instantaneous  and  the  time-average  measurements 
from  the  two  systems  during  periods  of  light  winds.  This  point  is  brought 
home  very  clearly  in  the  speed  and  direction  profiles  of  Figs.  50  and  51. 
When  wind  speeds  drop  below  2  m/s,  wind  direction  differences  become  large 
regardless  of  stability.  When  the  wind  speeds  are  larger,  the  agreement  be- 
tween the  rawinsonde  and  sodar  profiles  is  good.  The  two  cases  presented  here 
are  perhaps  more  spectacular  in  terms  of  the  wind  speed  effect  on  the  com- 
parison than  most  of  the  other  cases  examined.  Over  more  complicated  terrain, 
differences  in  speed  and  directions  between  sodars  and  rawinsondes  could  be 
much  larger.  Caution  must  be  exercised,  therefore,  in  interpreting  data  that 
compare  sodars  with  rawinsondes.  Nonetheless,  this  evaluation  does  indicate 
that,  under  proper  conditions,  reasonable  agreement  can  be  expected  between 
the  two  sets  of  measurements  as  though  both  techniques  measured  bulk  proper- 
ties of  the  wind  with  reasonable  accuracy. 
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8.  CHARACTERISTICS  OF  SODAR  w  SPECTRA 

The  spectra  of  the  vertical  wind  speed  derived  from  sodar  Doppler  measure- 
ments should,  in  principle,  correspond  to  spectra  from  the  sonic  anemometers, 
subject  to  the  effects  of  spatial  averaging  and  aliasing.  Spatial  averaging 
attenuates  fluctuations  with  scales  smaller  than  the  dimensions  of  the 
sampling  volume.  Aliasing  folds  energy  left  over  at  frequencies  above  n  ,  the 
Nyquist  frequency  (=  1/2  sampling  rate),  back  into  the  available  spectral 
bandwidth  (0  <  n  <  n  ).  In  the  presence  of  spatial  averaging,  the  energy 
folded  back  is  reduced  by  the  amount  lost  through  averaging.  A  schematic 
representation  of  the  distortions  introduced  on  a  typical  w  spectrum  for  two 
different  sampling  rates  is  given  in  Fig.  52(a).  In  this  example,  the  atten- 
uation from  spatial  averaging  is  assumed  to  commence  at  frequency  n,  =  0.02 
Hz.  The  wavelength  a,  corresponding  to  this  frequency  (a,  *  U/n. ,  where  U  is 
the  mean  horizontal  wind  component)  is  roughly  2tt  times  the  longest  dimension 
in  the  sampling  volume.  (A  sampling  volume  40  m  diam  x  40  m  long  is  assumed 
here  with  U  =  5  m/s.)  Because  of  the  sharp  spectral  attenuation  above  0.02 
Hz,  aliasing  is  confined  primarily  to  the  first  fold,  which  merely  raises  the 
energy  near  n  by  a  factor  of  2.  For  typical  beamwidths  used  in  most  sodar 
operations,  n  *  n,  *  0.02  Hz  at  100  <  z  <  300  m  for  moderate  wind  speeds. 

In  the  convective  boundary  layer  the  percentage  of  spectral  energy  con- 
tained in  frequencies  above  n  increases  as  height,  Z,  decreases.  Conse- 
quently the  uncertainties  in  the  observed  spectral  forms  and  in  the  measured 
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variances  also  increase  as  the  height  decreases.  Figure  52(b)  shows  the 
progression  of  the  spectrum  on  a  typical  day.  The  frequency  at  the  n  S  (n) 
spectral  maximum,  n  ,  is  nearly  constant  above  0.25  Z^ ,  (where  Z.  is  the 
boundary  layer  depth)  and  varies  inversely  with  height  below  that.  Within  the 
height  range  of  most  sodar  systems,  the  wavelength  at  the  spectral  peak  can  be 
approximated  by 


\n  =  U/nr 


6Z,     (Z  <  0.25  Z.) 
L 1.5Zi ,    (Z  >   0.25  Z.) 


(5) 


Spectral  energy  in  the  observed  bandwidth  also  drops  with  decreasing  Z. 

The  attendant  decrease  in  signal-to-noise  ratio  in  the  sodar  measurements 

serves  to  increase  further  the  uncertainty  in  the  spectral  and  variance  (a  ) 

estimates . 

In  the  stable  nocturnal  atmosphere,  the  w  spectral  scales  and  intensities 

are  more  strongly  controlled  by  stratification  than  by  Z.  Over  flat  terrain, 

within  the  stable  boundary  layer  (Kaimal  et  al.,  1972)  one  can  approximate 

x  using 
m    3 


-l 


^  .  Z(0.55  +  Z/L)  x  , 
«  L,  for  L  «  Z  , 


(6) 


where  L  is  the  Monin-Obukhov  length.  Within  the  height  range  of  our  com- 
parisons, \    would  be  roughly  an  order  of  magnitude  smaller  than  under 
unstable  conditions.  There  is  proportionally  less  energy  within  the  spectral 
bandwidth,  so  one  can  expect  to  find  larger  uncertainties  and  errors  in  the 
nighttime  spectra  than  in  the  daytime  spectra.  This  may  account  for  the 
increased  scatter  in  the  nighttime  <?w  values  in  Sec.  4.  An  improvement  in 
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Figure  52.  (a)  Schematic  representation  of  distortions  introduced  in 

the  w  spectrum  from  attenuation  due  to  spatial  averaging 
and  from  aliasing, 
(b)  Shift  in  spectral  behavior  with  height  and  its  implica- 
tions for  sampling  and  aliasing  errors. 
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accuracy  is  possible  in  the  presence  of  strong  gravity  waves  because  of  its 

large  contribution  to  variance  at  frequencies  below  n  . 

The  spectra  presented  in  Figs.  53  and  54  were  computed  from  time  series 

provided  by  AV.  No  significance  is  attached  to  the  choice  of  AV.  The  outputs 

are  treated  as  generic  signals  from  a  Doppler  sodar.  The  absence  of  liftup  at 

the  high  end  implies  extensive  influence  of  spatial  averaging  at  frequencies 

below  n  . 
o 

The  sonic  spectra  in  Fig.  54(a)  and  (b)  illustrate  the  effect  of  stability 
on  spectral  wavelengths  and  intensities  at  the  300  m  level.  The  sodar 
spectrum  shows  poor  agreement  with  the  sonic  spectrum  at  night;  spectral 
levels  are  greatly  enhanced.  The  high  a     levels  at  night  in  Sec.  4  can  now 
be  traced  back  to  this  distortion.  To  determine  how  much  of  this  distortion 
comes  from  aliasing,  the  sonic  time  series  was  converted  to  grab  samples  every 
24  s.  The  resulting  spectrum,  also  shown  in  Fig.  54(b),  has  the  same  shape  as 
the  sodar  spectrum,  but  one-half  the  energy. 

More  precise  estimates  of  the  contributions  from  aliasing  and  other  fac- 
tors, such  as  spatial  averaging  and  noise,  can  be  made  from  the  variances 
listed  in  Table  12.  Sonic  anemometer  variances  estimated  over  two  bandwidths, 
0-5  Hz  and  0-0.02  Hz,  are  listed  alongside  the  sodar  variances.  Sodar  vari- 
ances appear  to  be  10%  -  15%  lower  than  the  full  range  (5  Hz)  sonic  variances 
during  the  day  but  15%  -  20%  higher  than  the  sonic  variances  integrated  to 
0.02  Hz  (see  Table  12).  From  Table  13  (last  column)  we  find  the  sonic  anemom- 
eter variance  in  the  band  0.02  <  n  <  5  Hz  to  be  between  20%  and  25%  of  the 
total  variance  (0  <  n  <  5  Hz)  under  convective  conditions.  If  all  that 
variance  were  to  be  aliased  back  into  the  frequency  range  0  <  n  <  0.02  Hz,  the 

rati0  (aw\odar/(°w2)Son(0<n<5  Hz)  in  Tab1e  13  would  be  unity'  In  the 
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Figure  53.  Sodar  and  sonic  anemometer  w  spectra  at  (a)  150  m  and  (b) 
200  m  compared  for  morning  conditions. 
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Figure  54.  Sodar  and  sonic  anemometer  w  spectra  at  300  m  compared  for 
(a)  morning  and  (b)  nighttime  conditions. 
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Table  12.  Variances  of  w  from  the  sodar  and  the  sonic  anemometers 
for  the  cases  illustrated  in  Figs.  53  and  54 


Time      Date 

Height 

^°w  'sod 

(ow  >son 

(ow  Json 

(MST) 

(m) 

(0<n<0.02  Hz) 

(0<n<5  Hz) 

(0<n<0.02  Hz) 

0800-1000  9  Sept 

150 
200 
300 

0.846 
0.822 
0.762 

0.925 
0.951 
0.801 

0.711 
0.716 
0.644 

2220-2340  18  Sept 

300 

0.123 

0.082 

0.064 

Table  13.  Fraction  of  total  variance  sensed  by  the  sodar  and  sonic 
anemometers  over  bandwidth  0  <  n  <  0.02  Hz 


Time     Date      Height     (aw2V  .  (a,/), 


(MST)  (m) 


Jw   'sod  vow   ;son(0<n<0.02  Hz) 


/        2  \  o 

law   ;son(0<n<5  Hz)  Uw   )son(0<n<5  Hz) 


0800-1000  9  Sept  150  0.91  0.76 

0.75 


150 

0.91 

200 

0.86 

300 

0.95 

0.80 
2220-2340  18  Sept         300  1.50  0.78 
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unstable  case,  it  is  roughly  10%  smaller  at  all  three  heights,  which  suggests 
that  only  half  the  variance  above  0.02  Hz  is  aliased  back.  The  other  half 
represents  the  variance  lost  through  spatial  filtering.  In  the  stable  case  at 
300  m,  that  ratio  is  50%  higher  than  unity,  implying  contribution  from  sources 
other  than  the  true  w  fluctuations  in  that  spectral  bandwidth. 

The  time  series  used  for  computing  the  spectra  in  Fig.  54(b)  offer  clues 
as  to  the  cause  for  enhancement  in  the  sodar  spectrum.  Shown  from  top  to  bot- 
tom in  Fig.  55  are  the  fluctuations  in  w  recorded  by  the  sodar  (~0.04 
samples/s),  the  sonic  anemometer  sampling  slowed  down  to  the  same  rate  (~0.04 
samples/s),  and  the  original  sonic  anemometer  time  series  block-averaged  over 
10  s  nonoverlapping  blocks.  These  traces  are  counterparts  of  the  spectra  in 
the  0  <  n  <  0.02  Hz  range.  Not  surprisingly,  the  0.04/s-sampled  sonic  traces 
show  larger  excursions  than  the  block-averaged  traces.  In  the  sodar  traces, 
the  excursions  are  even  larger,  and  the  very  large  oscillation  that  occurred 
at  about  2230  MST  accounts  for  at  least  25%  of  the  excess  over  the  variance  in 
the  0.04/s-sampled  sonic  data.  Whether  this  oscillation  is  real  or  an  arti- 
fact of  data  processing  cannot  be  determined  now  since  the  only  other  w 
measurements  available  to  us  are  from  the  BA0  tower  600  m  away.  The  tower 
data  show  no  such  large  amplitude  oscillation.  At  this  time  one  can  only 
state  that  the  w  measurements  in  stable  air  have  to  be  edited  carefully  before 
any  a    value  is  accepted  as  correct. 
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Figure   55.     Time  series   of  w  corresponding  to  spectral    forms   shown   in 
Fig.   54(b)   in  the   frequency   range  0<n<0.02  Hz. 
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9.  CONCLUDING  REMARKS 

The  wind  measurements  analyzed  in  this  report  represent  the  state  of  the 

art  in  wind  sensing  with  commercially  available  Doppler  sodars.  Considering 

the  requirements  that  the  sodar  operation  be  unattended  except  for  maintenance 

and  repair,  and  that  the  data  be  subjected  to  no  editing  by  the  vendors,  the 

results  obtained  are  reassuring.  The  scatter  in  the  wind  speed  and  wind 

direction  data  compare  very  well  with  scatter  in  past  experiments  (Kaimal  et 

al.,  1980;  Gaynor  and  Korrell,  1981)  when  some  of  the  same  sodar  systems  were 

compared  at  the  BAO,  under  more  controlled  conditions.  In  these  data  some 

vendors  show  more  scatter  than  others,  but  much  of  that  can  be  attributed  to 

factors  such  as  rain,  high  winds,  cable  damage  (as  in  the  case  of  REM),  and 

transducer  failures  (as  with  XON).  If  data  from  these  suspected  periods  were 

eliminated,  the  different  systems  would  have  more  similar  scatter. 

The  measurement  of  a  with  sodars  seems  to  show  promise,  at  least  for 

w 

daytime  conditions.  Here  too,  vendor  performance  shows  variations  that  can  be 
attributed  in  some  cases  to  weather  and  equipment  failure,  but  individual  dif- 
ferences in  processing  the  data  might  also  be  a  factor.  The  predictable  beha- 
vior of  sodar  w  spectra  in  the  convective  boundary  layer  leads  us  to  believe 
that  a      measurements  can  be  made  to  within  10%  at  heights  above  100  m. 
However,  the  nighttime  results  are  not  so  encouraging.  More  work  is  needed  to 
ascertain  the  reasons  for  the  large  discrepancies  in  the  w  measurements  at 
night. 
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